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SUMMARY 

In  Part  I  of  this  study  (Lopez,  1988),  numerical  solutions  of  the 
axisymmetric  Navier-Stokes  equations  are  presented  and  compared  with 
results  from  experiments  for  a  confined  cylindrical  flow.  The  details  of  the 
vortex  breakdown  phenomenon  are  calculated  with  a  high  degree  of 
accuracy.  From  solutions  over  a  range  of  parameters  the  essential  features  of 
the  flow  are  obtained.  These  solutions  also  provide  flow  quantities  such  as  the 
vorticity  and  the  pressure  throughout  the  volume  which  would  be  difficult  to 
obtain  from  experiments.  In  this  paper  the  solutions  are  explored  and  the 
essential  physical  mechanisms  of  vortex  breakdown  in  this  particular 
geometry  are  identified.  These  mechanisms,  which  rely  on  the  production  of  a 
negative  azimuthal  component  of  vorticity  as  a  result  of  the  stretching  and 
tilting  of  the  predominantly  axially  directed  vorticity  vector,  are  elucidated 
with  the  aid  of  a  simple,  steady,  inviscid,  axisymmetric  equation  of  motion. 
This  equation  has  been  a  starting  point  for  most  studies  of  vortex  breakdown 
but  a  departure  in  the  present  study  is  that  it  is  explored  directly  and  not 
through  perturbations  of  an  initial  stream  function.  The  findings  are  then 
generalised  to  the  case  of  vortex  breakdown  in  swirling  pipe  flows. 
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1.  Introduction 

The  attraction  of  the  confined  swirling  flow  of  Part  I  to  a  study  of  vortex  breakdown 
is  that  the  flow  is  defined  by  only  four  variables,  H,  R,  ft  and  u,  each  of  which,  in  an 
experiment,  can  be  very  accurately  determined.  Similarly,  from  a  numerical  point  of 
view,  the  flow  is  confined  in  a  fixed  volume  with  very  well  defined  boundary  conditions. 
It  was  hoped  that  if  numerical  solutions  for  various  H/R  and  Re  =  QR2  jv  could  be 
obtained  which  very  accurately  predicted  the  resulting  flow  in  such  a  well  defined  ex¬ 
periment  then  the  physical  mechanisms  of  vortex  breakdown  could  be  elucidated  from 
these  solutions  and  further,  that  these  mechanisms  might  be  generalised  to  the  swirling 
flows  of  more  practical  significance  in  pipes  and  in  a  free  stream. 

Many  features  of  vortex  breakdown  have  been  well  recognised  and  the  early  works 
of  Squire  (1960)  and  Benjamin  (1962)  have  been  the  starting  point  for  many  subsequent 
studies.  A  departure  in  the  present  study  from  these  formulations  is  that  the  numerical 
solutions  are  explored  not  from  a  consideration  of  perturbations  to  an  initial  stream 
function  but  more  directly  from  the  development  of  the  azimuthal  component  of  vorticity. 
This  leads  to  criteria  which  depend  essentially  on  the  relationship  between  the  angle  of 
the  vorticity  vector  and  the  velocity  vector  on  stream  surfaces  upstream  of  breakdown. 

The  numerical  results  in  Part  I  demonstrate  the  accuracy  with  which  the  experi¬ 
mental  results  can  be  predicted  for  confined  flows  in  which  ‘vortex  breakdown’  develops 
and  this  paper,  Part  II,  explores  the  underlying  physical  mechanisms  which  account  for 
the  essential  features  of  the  flow  and  are  important  in  ‘vortex  breakdown’. 
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2.  Confined  swirling  flows 

2.1  Summary  of  the  principal  features  of  the  Sow. 

As  suggested  in  Part  I,  for  the  confined  flow  in  which  H/R  =  2.5  and  Re  >  1600, 
the  central  core  flow,  i.e  the  flow  returning  towards  the  rotating  endwall,  might  be  re¬ 
garded  as  practically  inviscid  apart  from  the  recirculation  zone  and  the  fluid  that  passes 
near  to  it.  This  first  approximation,  expected  to  be  increasingly  valid  as  Reynolds  num¬ 
ber  increases,  is  supported  by  Figure  1  which  shows  contours  of  the  streamfunction  V>, 
the  angular  momentum,  or  circulation,  T  and  the  total  head  H.  A  detailed  comparison 
between  these  contours  shows  that  in  this  central  core  region  T  and  H  are  approximately 
constant  on  stream  surfaces,  apart  from  those  surfaces  within  and  near  to  the  recircula¬ 
tion  bubble  where  F  and  H  are  relatively  small.  An  approximation  that  they  are  constant 
on  the  stream  surface  would  seem  unlikely  to  change  the  principal  features  of  the  flow. 
Of  course,  the  experiments  and  the  numerical  calculations  show  that  this  core  flow  is 
critically  dependent  on  Reynolds  number  but  mainly  in  the  sense,  it  seems,  that  Re 
determines  an  ‘upstream’  distribution  of  T(r)  and  H(r)  and  that  once  this  distribution 
is  established  the  subsequent  core  flow,  in  which  vortex  breakdown  occurs,  is  largely  an 
inviscid  rotational  flow  apart  from  the  flow  within  and  near  to  the  recirculation  zones. 

The  most  striking  features  of  this  core  flow  are  the  divergence  from  an  upstream 
narrow  core  to  a  much  larger  diameter  flow  which  is  almost  in  solid  body  rotation,  and 
correspondingly,  the  development  of  a  wave  as  a  result  of  an  ‘overshoot’  in  the  initial 
divergence.  The  apparent  effect  of  increasing  Reynolds  number  is  to  change  the  upstream 
distribution  of  T(r)  and  H(r)  in  a  way  such  as  to  increase  the  initial  divergence,  increase 
the  ‘overshoot’,  reduce  the  wavelength  of  the  resulting  wave  and  reduce  the  ‘damping’ 
of  this  wave.  The  contours  of  the  streamfunction  (Figure  2  from  Part  I)  suggest  that  the 
appearance  of  recirculation  bubbles  is  related  to  the  waviness  in  the  outer  flow.  These 
preliminary  observations  provide  useful  approximations  and  working  hypotheses  which 
are  explored  in  the  remainder  of  Section  2. 

2.2  Theoretical  considerations. 

Steady,  inviscid,  axisymmetric  swirling  flow  is  particularly  interesting  in  that  it  can 
be  looked  at  from  the  point  of  view  of  an  interaction  between  the  total  head  and  the 
angular  momentum  of  the  fluid  (both  of  which  are  conserved  on  a  stream  surface),  or,  in 
terms  of  a  balance  between  the  radial  pressure  gradient  and  the  centrifugal  force,  or,  in 
terms  of  the  generation  of  the  azimuthal  component  of  vorticity  through  the  stretching 
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and  tilting  of  vortex  lines. 

The  relationship  between  these  three  different  perspectives  is  as  follows.  In  terms 
of  the  radial  pressure  gradient  and  the  centrifugal  force,  t he  othei  radial  acceleration 
terms  are 

d u  du  i>2  I  dn 

(1) 


du  du  v2  1  dp 

UT  +  u’~a~  = - 

or  Oz  r  p  or 


In  terms  of  the  angular  momentum  (or  circulation)  T  —  rv ,  and  the  total  head  H  = 
p/p  +■  {u2  )-  v2  f-  w2)/ 2.  (1)  may  be  written  as 

du 


w 


dz 


Y_  0T 

r2  dr 


<9H  du- 

+  n-  — 

Or  Or 


(2) 


Thirdly,  since  the  azimuthal  component  of  vorticity  is  ?/  =  du/dz  -  dw/dr,  then  from 

(2) 

r  9r  du 

(3) 


r  ar  9H 

WT]  ~  r2  dr  ~  dr  ’ 


or,  since  the  streamfunction  ip  is  determined  by  u  —  - l/rdip/dz  and  w  —  l/rdip/dr 
and  since  T  and  H  are  constant  on  ipy  it  follows  from  (3)  that 


_  r  dr  (ffl 

^  r  dip  dip 


(4) 


Equation  (4)  has  been  a  starting  point  for  most  discussions  of  vortex  breakdown. 
Following  Squire  (1960),  it  is  usually  rewritten  by  replacing  tj  with  —  1  /r\/2ip.  Notwith¬ 
standing  that  T  and  H  are  functions  of  r  and  z,  the  equation  as  it  stands  is  of  such  a 
simple  form  that  we  have  pursued  its  implications  more  directly.  In  particular,  if  the 
curve  C  in  the  r,  z  plane  whose  radius  r  is  given  by  r  =  <r(z)  is  chosen  such  that  on  C 
the  streamfunction  is  a  constant,  i.e  ip(r,z)  =  ip\,  then  on  C 


or,  on  a  stream  surface 


and 


where  .4  =  r(V'i)r'(^i)  and  B  =  H'(^i). 


■■)  =  ^rVi)-<rH'r*i). 

( 7 

(5) 

dip  --  0  and  u  =  wcr'(z)-, 

(6) 

T](z)  =  -  -  5(7, 

(7) 

<7 

dp  ~  ~(^2  +  B^dcr, 

(8) 

=  H'(^,). 
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Equation  (7)  is  remarkable  since  it  provides  such  a  simple  expression  for  one  com¬ 
ponent  of  the  vorticity  in  terms  only  of  the  radius  of  the  stream  surface  yet  in  a  complex 
flow  where  the  vorticity  is  three  dimensional,  not  normal  to  the  velocity,  and  where  the 
stretching  of  vorticity  is  an  essential  mechanism  in  the  flow. 

Further  insight  into  the  implications  of  equations  (7)  and  (8)  is  gained  if  it  is 
assumed  that  at  some  upstream  station  zo  and  on  a  particular  stream  surface  il>  =  i/'i, 
the  radius  of  the  stream  surface  is  ao  and  the  azimuthal  and  axial  components  of  both 
the  velocity,  i.e.  r  and  w,  and  vorticity.  i.e.  r?  and  (.  are  known  and  have  the  values  vq, 
wq,  tj o  and  (o  (see  Figure  2).  In  this  case  A  and  B  may  be  determined  simply,  since 


r(  = 


dr  dT 

—  =  rw  — , 
dr  dil> 


so  that  on  i/> \ 


and 


r'(iM  =  — , 

u>0 


a  =  r(y>i)r'(vq)  = 

w0 


B  may  be  determined  from  (7)  evaluated  at  z  —  z0.  If  t/0  0,  then 

S  =  -i), 

<T0\wom  / 

and  if  7jo  =  0  then 


tfo  \w0tjo 
vo(o 


B  = 


<7q  Wq 


Thus  from  (7)  the  equation  for  tj  on  the  stream  surface  downstream  from  zq  is 

V  =  oc0J-° 

\  a  <7n/ 


for  tjo  =  0)  °r 

7/0  <7  \f3o)  <70  V/3o  )’  U 

for  7/o  y-  0,  where  ao  =  vo/wg  and  /3q  =  770/(0  are  the  tangents  of  the  helix  angle  for  the 

velocity  and  vorticity  respectively.  Thus,  downstream  of  zo,  the  ratio  of  7/  on  a  stream 
surface  to  its  value  on  the  surface  at  zo,  depends  only  on  the  ratio  of  the  tangents  of 
the  helix  angles  of  the  vorticity  and  velocity  on  the  surface  at  zo  and  on  the  ratio  of  the 
radius  of  the  stream  surface  to  the  radius  at  zo- 

Equation  (9)  has  interesting  implications.  Assuming  that  inviscid  behaviour  is  the 
dominant  mechanism  responsible  for  the  rapid  divergence  of  the  streamtubes  observed 
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in  vortex  breakdown,  a  reasonable  hypothesis  (in  the  absence  *~»f  bodies  in  the  flow  field) 
is  that  the  development  of  a  negative  azimuthal  component  of  vorticity  on  some  stream 
surfaces  is  necessary  if  the  axial  velocity  is  to  be  brought  to  zero  (see  Figure  2.)* 

If  this  is  so  then  from  equation  (9)  it  is  clear  that,  for  tj  q  positive,  rj  will  only  be'-ome 
negative  on  a  diverging  stream  surface  if  qo  >  /?p,  i  e.  vq/wo  >  r?o/Co-  A  helix  angle 
for  the  velocity  which  exceeds  that  of  the  vorticity  on  some  stream  surfaces  would  seem 
therefore  to  be  a  necessary  condition  for  vortex  breakdown  to  occur  since,  from  (9),  as 
a  I  cr  o  increases  from  unity,  then  for  op  fin.  decreases  to  zero  and  becomes  negative 

at  a  sufficiently  large  <t/ctq.  Figure  3  illustrates  the  behaviour  and  the  dependence  on 

ao/00- 

The  development  of  negative  tj  on  stream  surfaces  will  induce  a  negative  axial 
velocity  on  the  axis  which,  by  continuity,  will  lead  to  a  further  increase  in  a  and  corre¬ 
spondingly,  a  further  increase  in  negative  vorticity,  etc.  This  can  be  seen  as  a  kind  of 
‘positive  feedback’  which  accounts  for  the  relatively  rapid  divergence  of  stream  surfaces 
in  the  proximity  of  ‘breakdown’.  In  essence,  however,  this  is  a  matter  of  compatibility 
between  the  velocity  (U)  and  the  vorticity  (V  x  U)  arrived  at  from  a  consideration 
of  the  dynamics  of  vorticity  and  continuity.  It  follows  that  a  necessary  condition  for 
subsequent  vorioc  breakdown  (in  the  absence  of  viscous  or  turbulent  diffusion  and  the 
absence  of  bodies  in  the  flow  field)  would  seem  to  be  that  vq/u’o  >  770/Co 

A  reduction  in  azimuthal  vorticity  with  increasing  radius  is  readily  understood 
physically  for  the  simple  case  of  a  flow  which  is  cylindrical  upstream  and  has  tj  =  0.  If 
two  material  points  on  a  vortex  line  (which  is  axial  in  this  cylindrical  flow)  move  with 
the  fluid  and  the  leading  point  advances  to  a  location  where  the  radius  has  begun  to 
increase,  then  this  leading  point  will  experience  a  reduction  in  azimuthal  velocity  due 
to  the  conservation  of  T  on  the  stream  surface.  The  material  line  between  the  points  is 
therefore  stretched  and  tilted,  and  since  the  vorticity  moves  with  this  material  line,  a 
negative  azimuthal  component  of  vorticity  is  generated.  This  accords  of  course  with  the 

*  For  simplicity,  in  this  paper  z  is  taken  to  increase  in  the  direction  of  the  primary 
axial  velocity  in  the  core  flow  and  correspondingly  negative  azimuthal  vorticity  will 
induce  an  axial  velocity  in  the  opposite  direction.  In  Part  I  increasing  z  is  in  the 
opposite  direction  away  from  the  rotating  end  wall  and  the  sign  of  the  corresponding 
azimuthal  vorticity  is  the  reverse  of  its  sign  in  this  paper. 
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tjo  =  0  case  of  equation  (9).  More  generally  however,  for  tj0  y  0,  tlien  from  (8)  or  (9), 

dv  _  zUlf  (^1  ,  ,\£0 

drr  0-o  \  \  cr2  )  Bo 

which  implies  not  a  reduction  but  an  increase  in  azimuthal  vorticity  with  increasing 
radius,  i.e,  in  the  absence  of  a  body  in  the  flow  field,  ‘negative  feedback’  in  terms  of  the 
above  compatibility  argument,  f 

(% 

\  ex¬ 
it  follows  that  for  (x/cxq  >  1  and  rj0  positive  then  for  cto/Bo  <  0.5,  dr]/d.a  >  0;  for 
0.5  <  a o//3q  <  1,  X]  remains  positive  and  at  sufficiently  large  cr j a o,  d-qjdcr  >  0;  and 
only  for  ao/Bo  >  1,  does  p  become  negative  for  sufficiently  large  ex/ao-  (See  fig.  3). 
[This  more  general  result  is  not  easily  interpreted  physically  but  from  an  extension  of 
the  above  simple  argument  for  the  case  of  q  =  0,  one  can  see  that  if  Bo  >>  cto  it  is 
plausible  that  the  vorticity  vector  could  be  rotated  to  an  even  larger  helix  angle  by  a 
change  in  the  velocity  vector  with  radius  (u  and  w  decrease  with  increasing  r)  so  that 
the  result  that  dx]jd<x  >  0  for  Bo  >  2a<j  is  not  as  surprising  as  it  may  seem  at  first  sight.) 

An  analogy  in  these  terms  with  supersonic  flows  seems  possible  in  the  sense  that 
strong  negative  feedback  (drj/da  positive)  suggests  stability  and  an  absence  of  waves  on 
the  stream  surface.  There  is  not  a  simple  connection,  however,  between  this  result  and 
the  supercritical,  subcritical  distinction  described  by  Benjamin  (1962,  1967). 

With  these  general  considerations  in  mind,  the  vortex  breakdown  features  of  the 
confined  swirling  flow  may  be  interpreted  as  follows. 

2.3  Physical  mechanisms. 

The  summary  in  Section  2.1  and  these  theoretical  considerations  suggest  that  an 
important  question  to  ask  is  why  the  strong  vortical  core  flow  begins  to  diverge?  A 
simple,  if  simplistic  answer  is  clear  from  equation  (1).  Downstream  from  the  point  where 
u  =  0  and  the  stream  surface  has  its  smallest  radius  (or,  more  generally,  downstream 
from  a  region  in  which  the  flow  is  cylindrical)  there  will  be  divergence,  i.e.  positive 
duj dz  only  if  the  centrifugal  force  exceeds  the  radial  pressure  gradient  [in  the  absence  of 
viscous  stresses).  For  the  confined  flow,  this  imbalance  clearly  arises  from  the  fact  that 
the  radial  distribution  of  T  and  H  in  the  narrow  upstream  core  region  is  determined  by 
the  ‘upstream’  history  of  a  fluid  particle  as  it  travels  on  its  closed  stream  surface.  In 


particular,  the  particle  acquires  and  loses  most  of  its  angular  momentum  and  total  head 
through  the  action  of  viscous  stresses  in  the  boundary  layers  near  the  surfaces  of  the 
volume  (Figure  1).  There  is  no  reason  why  the  radial  distributions  of  T  and  H  should 
lead  to  t]  -  0  (equation  (3)),  or  to  cvclostrophic  balance  (equation  (1)).  Quite  the 
contrary,  in  fact,  since  the  turning  of  the  flow  from  radially  inward  towards  the  lower 
rotating  endwall  requires  the  centrifugal  force  to  exceed  the  radial  pressure  gradient. 
This  is  supported  by  Figure  4  which  shows  the  calculated  contours  of  t r/r  -  l/pdp/dr 
for  the  case  Re  -  1994.  H  R  -  2.5.  (These  contours  provide  an  interesting  view  of  the 
force  field;  they  are  not  surprising  in  the  sense  that  they  accord  with  the  waviness  of  the 
stream  surface  and  the  corresponding  sinusoidal  like  variation  in  the  radial  acceleration.) 

The  reason  for  the  initial  ‘overshoot’  in  the  divergence  of  the  stream  surfaces  and 
the  subsequent  waviness  can  be  seen  as  follows.  From  (4),  we  have 

du  TdT  dH  dw 

^  =  +  (10) 
For  the  confined  swirling  flow,  If  z  =  20  where  u  —  0  and  du/dz  has  a  local  maximum, 
i.e.  where  a  stream  surface  has  its  smallest  radius,  then  on  a  particular  stream  surface 
of  radius  r  =  <r{z),  whose  radius  is  cr  =  2  =  zq, 


where  .4  and  B  are  the  particular  values  of  FdT /dx{;  and  dH/d if;  on  this  stream  surface. 

The  characteristic  shape  of  the  stream  surface  downstream  from  20  can  be  explored 
as  follows.  Since  u  is  wcr'  (equation  (6)),  equation  (11)  may  be  approximated  by 

w{<t)<t"  -  w(<ro)cr'' {zq)  ^  a(  1  -  1  -  B(cr  -  cr0)  +  -  .  (12) 

\<7  (Xq  /  UT  (j  UT  a ,1 

(Note  that  ui(cr)  =  ui(o-o)  corresponds  with  a  small  perturbation  linearization).  Immedi¬ 
ately  downstream  of  20  for  example,  the  diverging  stream  surface  implies  an  increasing  a 
and  cr"  is  expected  to  decrease,  i  e.  cr"(z)  <  cr'^co)  because  .4(  1/cr -  l/Vo )  and  -F(cr  — ctq) 
are  both  negative,  assuming,  in  accord  with  Section  2.2,  that  since  breakdown  occurs, 
ao  >  0o  and  B  is,  therefore,  positive.  Note  however  that  dw/dr\a  -  dw/dr\(Jn  is  expected 
to  be  positive  because  diverging  stream  surfaces  imply  by  continuity,  a  reduction  in  w 
and  a  relatively  larger  reduction  in  w  is  expected  near  r  =  0  than  at  large  r,  (dw/dr |ffn 
is  expected  to  be  negative).  Thus,  downstream  of  20,  cr"  is  expected  to  decrease,  but  the 
rate  at  which  this  occurs  will  depend  on  w.  The  continued  divergence  of  the  flow  and 
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corresponding  reduction  in  a"  tnerefore  leads  to  a  radius  of  the  stream  surface  at  which 
cr"  =  0  (infinite  radius  of  curvature  of  the  stream  surface).  At  this  point  the  local  slope 
<r\  is  a  maximum.  Since  the  flow  continues  to  diverge  downstream  from  this  point,  cr" 
becomes  negative  and  the  stream  surface  begins  to  reduce  in  slope,  turning  back  towards 
the  axis. 

An  alternative  explanation  for  the  flow  divergence  and  subsequent  turning  of  the 
flow  towards  the  axis  may  be  sought  in  terms  of  equation  (1).  An  imbalance  in  v2/r  — 
2/ pdpjdr  gives  rise  to  a  positive  value  of  cr" .  On  this  diverging  stream  surface,  v'/r  is 
decreased  due  to  the  conservation  of  T  and  at  the  same  time,  w  and  v  are  reduced  due  to 
the  divergence,  leading  to  an  increased  pressure  with  relatively  larger  changes  for  small 
r  and  a  reduction  in  the  radial  pressure  gradient.  While  not  self-evident,  the  reduction 
in  v2/r  must  exceed  the  reduction  in  dp/ dr  since,  as  outlined  above,  tr"  approaches  zero 

A  further  alternate  and  more  direct  explanation  may  be  sought  in  terms  of  the  kine¬ 
matics  of  vorticity  (equation  (4)).  As  discussed  in  Section  2.2,  for  positive  B,  (cto  >  0o) 
the  diverging  flow  gives  rise  to  a  reduction  of  the  positive  azimuthal  component  of  vor¬ 
ticity  and  the  development,  with  increasing  radius,  of  negative  azimuthal  vorticity.  In 
this  case,  the  ‘positive  feedback’  discussed,  leads  to  the  further  development  of  nega¬ 
tive  azimuthal  vorticity  and  the  divergence  is  expected  to  continue  until  the  negative 
azimuthal  vorticity  grows  to  a  sufficiently  large  magnitude  for  it  to  turn  the  diverging 
flow  back  towards  the  axis. 

The  fact  that  the  divergence  ‘overshoots’  and  that  a  wave  in  the  stream  surface 
results  may  be  approximately  seen  by  linearising  equation  (12)  about  the  radius  a  =  cr\ 
at  which  a"  =  0.  Thus  for  cr  —  cr i  +  £(z),  then  from  (12) 

~  (13) 

where  perturbations  in  dwfdr  have  been  neglected  compared  with  du/dz.  A  wave  in  the 
stream  surface  is  therefore  established  with  wave  number  k  given  by  w\k2  ~  A/ir\  +  B. 

The  shape  of  a  characteristic  stream  surface  in  the  core  flow  downstream  of  zo  can 
be  summarized  as  follows.  The  stream  surface  c t(z)  diverges  from  <ro  to  crj  and  cr"  is 
reduced  from  a  maximum  to  zero.  At  the  radius  a  =  oq,  the  positive  slope  cr'  is  a 
maximum  and  the  flow  continues  to  diverge,  £  >  0,  but  a" ,  or  ,  becomes  negative  and 
the  surface  begins  to  turn  towards  the  axis.  This  overshoot  of  the  radial  displacement 
leads  to  a  subsequent  oscillation  in  the  displacement  of  the  stream  surface  about  cr  —  cr- 
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i.e.  (  —  0.  The  wavelength  of  this  displacement  is  approximately 


A 


~  2ww^'5 


(T_dT 
\«r^  dtp 


dH\-°5 

dtp) 


(14) 


where  w\  is  the  axial  velocity  at  r  =  oq. 

Using  the  expressions  for  A  and  B  from  Section  2.2,  the  wavelength  can  also  be 
approximated  by 


(15) 


A  comparison  between  the  wavelength  evaluated  from  (15)  and  that  determined  from  the 
periodicity  in  the  numerical  solution  is  presented  in  Table  1.  The  wavelength  predicted 
by  equation  (15)  was  evaluated  on  a  stream  surface  where  T  and  H  are  approximately 
constant  (selected  by  comparing  the  contours  of  T,  H  and  ip  as  in  Figure  1).  The 
wavelength  from  the  numerical  solution  was  determined  from  the  contours  of  v2/r  — 
l/pdp/dr.  Of  course,  the  wavelength  found  in  both  cases  can  be  regarded  only  as 
a  typical  wavelength.  As  expected  the  agreement  (approximate)  in  Table  1  between 
the  inviscid  prediction  of  equation  (15)  (and  equation  (14))  and  the  numerical  solution 
increases  with  increasing  Reynolds  number.  The  difference  between  the  predictions  of 
equations  (14)  and  (15)  reflects  the  fact  that  equation  (14)  is  evaluated  in  the  region  of 
the  waves,  where  the  flow  is  most  inviscid,  whereas  equation  (15)  assumes  inviscid  flow 
from  zo  where  the  radius  of  the  stream  surface  is  a  minimum.  Closer  examination  of 
Figure  1  shows  that  the  approximation  that  T  and  H  are  constant  on  stream  surfaces  is 
not  accurate  in  this  legion.  Particularly  at  lower  Reynolds  numbers  it  is  not  surprising 
therefore  that  the  azimuthal  viscous  stresses  should  lead  to  some  additional  reduction  in 
the  azimuthal  component  of  vorticity  in  this  region  and  correspondingly  that  equation 
(9)  would  predict  a  larger  radius  than  the  numerical  solution,  for  a  given  reduction  in  77. 
This  is  consistent  with  the  prediction  from  equation  (15)  of  a  shorter  wavelength  than 
is  four.d  from  the  numerical  solution. 

As  a  check  whether  the  model  describes  the  principal  features  of  the  flow,  how¬ 
ever,  the  agrt.-  r. ■  seems  satisfactory.  In  accord  witli  this  discussion  and  that  of  Section 
2.2  the  subsequent  development  of  a  second  recirculation  zone  could  be  expected  if  the 
wave  is  of  s  'fficit-w  amplitude  for  the  subsequent  diverging  stream  surface  to  develop 
sufficient  negative  azimuthal  vorticity  to  again  bring  the  axial  velocity  to  zero.  The 
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converging  stream  surface,  by  contrast,  generates  positive  azimuthal  vorticity  which 
accelerates  the  axial  iiow. 

The  principal  features  and  physical  mechanisms  of  the  ‘breakdown’  in  the  confined 
swirling  flow  are  therefore  comprehensible  in  terms  of  inviscid  phenomena.  The  initial 
divergence  from  a  narrow  core  is  the  result  of  the  upstream  radial  distribution  of  H  and 
T,  the  subsequent  waviness  of  the  stream  surface  is  a  result  of  a  stationary  inertial  or 
centrifugal  wave  and  the  development  of  one  or  more  recirculation  zones  is  a  consequence 
of  the  generation  of  negative  azimuthal  vorticity  through  the  stretching  and  tilting  of 
vortex  lines  in  a  diverging,  swirling  inviscid  flow.  At  a  Reynolds  number  above  that 
at  which  breakdown  first  occurs,  the  apparent  ‘critical’  dependence  of  the  flow  on  Re, 
discussed  in  Part  I,  arises  essentially  from  the  dependence  of  the  radial  distributions 
of  T  and  H  (in  the  core  flow)  on  the  action  of  viscous  stresses  in  the  boundary  layer 
regions  near  the  walls.  At  higher  Reynolds  number  the  core  flow  behaves  in  a  broadly 
inviscid  manner  and  the  recirculation  zones  appear  to  have  relatively  little  dynamical 
significance,  except  insofar  as  they  are  within  the  region  of  strong  negative  azimuthal 
vorticity. 

These  mechanisms  which  appear  to  account  for  the  principal  features  of  breakdown 
in  the  confined  flow  led  to  our  considering  their  applicability  to  swirling  flows  in  pipes. 
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3.  Swirling  flows  in  pipes  and  in  a  free  stream 

3.1  Preliminary  considerations. 

The  relatively  sudden  appearance  of  a  rapid  divergence  in  the  stream  surfaces,  and 
the  corresponding  appearance  of  recirculation  zones,  occurs  in  some  swirling  pipe  flows 
and  in  the  vortex  breakdown  observed  in  flows  over  delta  wings  at  high  angle  of  attack. 
A  comparison  with  the  ‘breakdown’  observed  in  the  confined  swirling  flow  suggests  that 
these  breakdown  regions  are  qualitatively  similar,  at  least  up  to  the  first  recirculation 
region.  An  essential  difference  however,  is  that  for  swirling  pipe  flows  and  to  some  extent 
the  delta  wing  flow,  the  ‘upstream’  region  may  be  regarded  as  essentially  a  cylindrical 
flow  for  which  initially  v2 /r  =  l/pdp/dr  or,  in  terms  of  equation  (4) 


du  TdT  dH  dw 
dz  r  dip  dip  +  dr 

For  this  case,  the  question  that  is  immediately  posed  is  how  does  ‘breakdown’  begin? 
One  of  us  (Lopez)  has  obtained  time  dependent  numerical  solutions  of  the  Navier  Stokes 
equations  for  various  swirling  flows  in  a  pipe  by  the  methods  outlined  in  Part  I.  The 
particular  family  of  flows  that  has  been  calculated  is  given  at  t  =  0  by  a  zero  radial 
component  of  velocity,  an  azimuthal  component  of  velocity 


v  = _ —1 _ (i_e-r’/°  ) 

rU-e-V1  >' 

and  an  axial  component  of  velocity 


(18) 


v>  =  W„(l  +  Wet-’’!''). 


(17) 


These  distributions  have  been  found  to  be  good  fits  to  various  experimental  observations 
of  vortices  undergoing  breakdown  (eg.  Garg  &  Leibovich,  1979).  This  azimuthal  velocity 
distribution  has  the  further  merit  that  it  has  the  same  form  as  the  exact  solution  for  the 
diffusion  of  vorticity  from  a  line  vortex  (Burgers,  1940). 

The  distributions  of  velocity  in  equations  (16)  and  (17)  have  the  corresponding 
azimuthal  and  axial  components  of  vorticity, 

2WoaWcr  _riui 
T'=— T— e  '  - 


c  = 


2v<= 

a(l  —  e”1) 


and 
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and  in  terms  of  the  helix  angles  for  velocity  and  vorticity  discussed  in  Section  2.2,  on 
stream  surfaces  at  z  =  0 

=  /  g  \2Vc3*2  (l-e-r^) 

0o  Ke-ll  iycr2(l  +  ^Vce-rl/«,)’ 

where  14  =  vc/We c.  Within  the  range  1  <  Wc  <  2,  a/0  does  not  vary  substantially  for 
stream  surfaces  for  which  0  <  r/a  <  1  and  at  r/a  =  1 

oo  =  Wce2  V? 

00  ~  (e  -  1  )(Wc  +  e)W}' 

This  relationship  between  a<j  and  0q  a.t  r/a  =  1  for  various  values  of  Ve  and  Wc  is  shown 
in  figure  5  and  the  regime  is  marked  for  which  the  inviscid  criterion  for  vortex  breakdown 
(ao  >  0a)  is  met.  The  points  where  numerical  solutions  have  been  obtained  are  also 
given  in  the  figure.  For  a  particular  flow  to  be  discussed,  14  =  1.5,  We  =  1.25  and  at 
r  =  a,  ao/0o  is  1.95. 

The  initial  conditions  for  the  calculated  flows  are  that  equations  (16)  and  (17)  are 
the  velocities  at  t  =  0  for  all  z.  The  upstream  boundary  condition  is  that  these  velocities 
are  imposed  for  t  >  0.  The  radial  boundary  condition  at  r  =  R  is  u  =  0,  dv/dr  =  0  and 
t]  =  0.  The  usual  outflow  condition,  d/dz  =  0  is  applied.  (A  further  discussion  of  these 
boundary  conditions  and  their  implementation  is  set  out  in  Appendix  I.)  With  these 
initial  and  boundary  conditions  the  presence  of  viscosity  is  essential  for  the  evolution  of 
this  flow  -  inviscid  flow  would  not  change  with  time.  An  example  of  the  development 
of  the  flow  for  a  particular  Reynolds  number  is  shown  in  Figure  6.  Of  course,  at  very 
low  Reynolds  number,  viscous  diffusion  dominates  the  flow  development.  At  higher 
Reynolds  numbers  however,  we  have  found  that  for  the  same  initial  cylindrical  velocity 
distributions,  the  time  taken  before  a  recirculation  zone  appears  increases  with  Reynolds 
number  (defined  as  Re  =  W^af  u).  Unlike  the  confined  swirling  flow,  for  these  flows 
it  is  the  diffusion  of  axial  vorticity  which  leads  initially  to  u2/r  —  1/pdp/dr  becoming 
positive.  Alternatively,  we  anticipate  that  it  must  be  the  diffusion  of  vorticity  which 
leads  to  a  radial  redistribution  of  T  and  to  the  stretching  and  tilting  of  vortex  lines  (due 
to  the  axial  change  in  with  a  corresponding  reduction  in  the  initial  positive  azimuthal 
component  of  vorticity  with  axial  distance,  and  the  subsequent  beginning  of  an  ‘invisdd’ 
breakdown  process. 

This  broad  view  of  the  likely  role  of  viscosity  and  the  development  of  azimuthal 
vorticity  is  illustrated  in  the  calculations  of  v,  rf  and  i/>  shown  in  Figure  6  and  leads  to 
the  following  more  detailed  considerations. 
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3.2  Physical  Mechanisms. 

In  considering  the  results  shown  in  Figure  6,  particular  attention  is  drawn  to  the 
reduction  of  the  azimuthal  velocity  and  vorticity  with  distance  downstream.  Initially, 
this  reduction  is  due  to  viscous  effects,  but  it  is  clear  that  by  t  =  40,  there  is  a  slight 
divergence  of  streamlines  which,  by  the  inviscid  mechanisms  of  section  2,  contributes  to 
this  reduction.  At  t  =  80  this  azimuthal  vorticity  has  been  reduced  to  zero  on  stream 
surfaces  at  a  downstream  distance  of  about  20a  and  out  to  a  radius  of  about  a  and  a  fur¬ 
ther  divergence  of  these  stream  surfaces  has  generated  a  negative  azimuthal  component. 
The  positive  feedback  mechanisms  described  in  section  2.2  drive  this  further  develop¬ 
ment  of  a  negative  azimuthal  component  leading,  by  t  =  240,  to  a  small  recirculation 
zone  on  the  axis,  rapid  changes  in  the  azimuthal  vorticity  ahead  of  this  region  where 
the  streamlines  diverge,  and  the  evident  propagation  upstream  of  the  region  of  negative 
azimuthal  vorticity  due  to  its  own  induced  velocity.  The  subsequent  appearance  of  a 
second  recirculation  zone  located  at  a  distance  downstream  where  the  outer  stream  lines 
diverge  and  where  the  azimuthal  component  of  the  vorticity  has  a  second  maximum  is 
consistent  with  the  establishment  of  a  wave  on  the  outer  stream  surfaces  by  essentially 
the  same  mechanism  as  in  the  confined  swirling  flows. 

Further  insight  into  the  mechanics  of  the  flow  is  obtained  from  a  detailed  comparison 
between  the  results  shown  in  figures  7(b),  7(c),  and  7(d).  Figures  7(a)  and  7(b)  show 
the  development  of  the  flow  for  which  Vc  =  1.75,  Wc  =  1.6  and  Re  =  300  (with  a 
corresponding  cto/A)  =  1.91)  from  t  =  227  to  t  =  250.  Figure  7(c)  shows  a  corresponding 
development  in  an  unphysical  case  in  which,  for  the  above  flow,  at  t  =  227  the  viscosity 
is  suddenly  doubled,  (i.e.  Reynolds  number  reduced  from  300  to  150).  Figure  7(d)  is 
the  case  in  which  the  viscosity  is  suddenly  halved  (Re  =  600)  at  t  =  227.  A  comparison 
between  figures  7(b)  and  7(c)  shows  that  the  subsequent  effect  of  a  sudden  increase 
in  the  viscosity  is  to  diffuse  the  axial  vorticity  and  increase  the  initial  divergence  of 
streamlines  but  reduce  the  magnitude  of  the  maximum  negative  component  of  azimuthal 
vorticity  from  -1.84  to  -1.37,  to  reduce  the  size  of  the  recirculation  bubble,  and  reduce 
the  rate  at  which  the  field  of  negative  vorticity  (and  the  bubble)  propagate  upstream.  By 
contrast,  comparison  between  figures  7(b)  and  7(d)  shows  that  reducing  viscosity  reduces 
the  diffusion  of  axial  vorticity  which  reduces  the  initial  divergence  of  streamlines;  the 
presence  of  the  bubble  and  divergence  of  the  stream-lines  immediately  ahead  of  it  and  the 
reduced  diffusion  of  axial  vorticity  increases  the  magnitude  of  the  maximum  negative 
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component  of  azimuthal  vorticity  that  is  reached  from  -1.84  to  -2.17.  The  bubble  is 
correspondingly  larger,  it  propagates  upstream  a  little  more  quickly,  the  amplitude  of 
the  wave  on  the  outer  stream  surface  is  greater  and  a  small  second  recirculation  zone 
begins  to  emerge. 

These  results  are  quite  consistent  with  the  discussion  in  section  2.2.  For  this  flow 
the  necessary  criterion  for  breakdown  to  occur  (ao  >  0q)  is  satisfied.  The  development 
of  the  flow  requires  viscosity  to  initiate  a  reduction  of  the  azimuthal  components  of 
vorticity  and  velocity  with  distance  downstream  and  to  initiate  the  divergence  of  the 
streamlines.  However,  the  development  of  the  breakdown  is  essentially  an  invisdd  pro¬ 
cess  driven  by  the  positive  feedback  mechanisms  described  in  section  2.2.  The  reduction 
in  viscosity  introduced  after  a  bubble  has  appeared,  figure  7(d),  shows  the  dominant  role 
that  this  invisdd  mechanism  plays.  In  some  respects  the  phenomenon  in  this  swirling 
pipe  flow  case  is  not  unlike  boundary  layer  transition  in  the  general  sense  that  the  ini¬ 
tial  development  (Tollmein-Schlichting  waves  in  boundary  layer  transition)  depends  on 
viscosity  and  the  final  transition  process  is  a  consequence  of  ‘positive  feedback’  in  the 
generation,  in  that  case,  of  a  longitudinal  component  of  vortidty. 

3.3  Comparison  with  experiment. 

While  more  detailed  comparisons  with  experiment  over  a  wide  range  of  parameters 
have  not  yet  been  made,  the  essential  features  of  the  evolution  of  an  axisymmetric  vortex 
breakdown  are  evident.  The  description  by  Sarpkaya  ( 1971)  of  the  initial  axisymmetric 
swelling  of  the  vortex  core  together  with  Escudier’s  (1986)  evolutionary  sequence  in  his 
Figure  9.12  (reproduced  here  in  Figure  8)  are  seen  in  the  calculation  of  the  axisymmetric 
evolution  of  Figure  6.  Noting  that  in  the  calculations  the  time  scale  of  the  evolution 
is  given  by  Woo /a,  and  recalling  that  the  Reynolds  number  is  based  on  the  core  radius 
a  rather  than  the  pipe  diameter,  we  find  that  the  time  scale  of  the  evolution  compares 
favourably,  (at  a  comparable  Re)  with  that  obtained  in  Escudier’s  (1986)  experiment. 
(Due  to  the  lack  of  a  quantitative  measurement  in  the  experiment,  only  a  rough  esti¬ 
mate  of  the  core  diameter  a  could  be  made).  A  comparison  between  the  structure  of  the 
leading  bubble  with  the  flow  visualizations  of  Escudier  (1986)  shows  good  agreement. 
For  comparable  times  in  both  the  experiment  and  the  calculation,  the  bubble  migrates 
upstream  approximately  the  same  distance  (four  bubble  diameters)  from  the  point  of 
initial  swelling  of  the  core.  The  appearance  of  a  second  breakdown  zone  at  about  one 
bubble  diameter  downstream  from  the  leading  bubble  is  evident  in  both  experiment  and 
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calculation.  In  the  calculation  this  second  recirculation  zone  is  very  unsteady  by  com¬ 
parison  with  the  leading  bubble  and  the  flow  visualization  is  suggestive  of  asymmetries 
in  this  region  of  the  flow,  whereas  the  leading  bubble  is  evidently  largely  axisymmetric. 
Also  very  evident  in  the  computed  leading  bubble  is  the  ‘two  celled’  structure  first  iden¬ 
tified  in  bubble  breakdowns  by  Faler  k  Leibovich  (1978).  The  main  outer  cell,  or  vortex 
ring,  has  an  associated  negative  azimuthal  component  of  vorticity  and  the  weaker  inner 
ring  has  a  positive  azimuthal  component  of  vorticity. 
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4.  Concluding  remarks 

Many  of  the  features  of  vortex  breakdown  have  been  well  recognised  previously. 
The  early  works  of  Benjamin  (1962,  1967)  have  been  the  starting  point  for  many  sub¬ 
sequent  studies.  However,  for  the  confined  flow  and  the  swirling  pipe  flow,  we  have  not 
found  the  ‘conjugate  flow’  theory  of  particular  assistance  in  illuminating  the  mechanics 
of  these  flows,  for  the  following  reasons.  The  theory  assumes  a  cylindrical  flow  upstream, 
whereas  in  the  confined  flow,  there  is  no  region  far  upstream  of  breakdown  in  which  the 
centrifugal  forces  balance  the  radial  pressure  gradient.  Similarly,  there  is  no  evidence  of 
the  importance  of  energy  dissipation  in  the  outer  flow.  In  the  swirling  pipe  flow  case, 
whereas  the  ‘conjugate  flow’  theory  assumes  an  inviscid  development,  the  calculations 
and  physical  mechanisms  for  the  generation  of  negative  azimuthal  vorticity  that  we  have 
considered  depend  initially  on  the  diffusion  of  axial  vorticity.  It  seems  nevertheless  true 
that  having  established  breakdown,  the  azimuthal  vorticity  which  is  generated  propa¬ 
gates  upstream  on  an  inertial  time  scale  and  the  diffusion  of  vorticity,  which  is  initially 
responsible  for  the  breakdown  is  of  lesser  significance  and  continues  to  operate  on  the 
much  slower  viscous  time  scale. 

From  the  discussion  in  this  paper,  the  essential  feature  of  vortex  breakdown  is  the 
generation  of  negative  azimuthal  vorticity.  For  the  confined  flow,  the  rapid  generation 
of  this  component  of  vorticity  is  a  consequence  of  the  upstream  distribution  of  T(r)  and 
H(r)  and  the  non- cylindrical  nature  of  the  upstream  flow.  For  the  swirling  pipe  flow,  we 
find  the  initial  reduction  in  the  azimuthal  component  of  vorticity  to  be  due  to  viscous 
diffusion  and  the  resultant  stretching  and  tilting  of  axial  vorticity.  As  the  flow  begins 
to  diverge  however,  the  further  production  of  nega  ive  azimuthal  vorticity  is  dominated 
by  inviscid  mechanisms  and  there  is  a  close  similarity  between  the  confined  flow  and  the 
swirling  pipe  flow  in  this  region. 

A  preliminary  application  of  these  ideas  to  the  more  complex  problem  of  vortex 
breakdown  in  delta  wing  flows  suggests  that  the  radial  and  axial  velocity  distributions 
in  the  core  flow  which  establish  a  characteristic  ao/^0  for  each  angle  of  attack  will  play 
a  critical  role  in  determining  whether  breakdown  can  occur  and  if  so  the  strength  and 
location  of  it.  Similarly,  by  broad  analogy  with  the  swirling  pipe  flow,  turbulent  diffusion 
could  be  expected  to  be  important  if  the  flow  initially  is  in  approximate  cyclostrophic 
balance.  In  this  case  increased  turbulent  diffusion  would  be  expected  to  reduce  the  size 
of  the  breakdown  bubble  and  to  reduce  the  distance  to  breakdown.  Of  course,  an  outer 
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flow  which  creates  diverging  stream  surfaces  will  shorten  the  distance  to  breakdown 
and  increase  the  relative  importance  of  the  inviscid  mechanism  in  reducing  an  initially 
positive  helix  angle  of  the  vorticity.  [This  effect  of  a  diverging  outer  flow  has  been 
observed  in  swirling  pipe  flows  and  may  be  seen  in  a  comparison  between  the  diverging 
pipe  results  of  Sarpkaya  (1971)  and  the  constant  diameter  results  of  Harvey  (1962)]. 

The  essential  ideas  which  link  the  generation  of  negative  azimuthal  vorticity  with  the 
helix  angles  of  the  velocity  and  vorticity  vectors  will  be  applicable  also  to  swirling  flows 
past  axisymmetric  bodies  located  on  the  swirl  axis.  The  positive  feedback  mechanism 
could  account,  for  example,  for  the  sensitivity  of  these  flows  to  the  intrusion  of  external 
probes  leading  in  some  cases  to  a  rapid  development  of  vortex  breakdown.  For  ao  >  0o 
the  divergence  of  stream  surfaces  over  a  body  could  be  expected  to  lead  to  a  forward 
recirculation  zone  ahead  of  the  body  if  ao//?o  is  sufficiently  large  or  the  body  has  a 
large  enough  diameter.  Such  a  flow  with  a  forward  recirculation  zone  has  been  observed 
recently  by  Joubert  (private  communication).  As  discussed  in  section  2.2,  the  flow  for 
the  case  where  (3q  >  2ao  is  much  less  clear;  the  negative  feedback  that  would  accompany 
an  increase  in  radius  of  a  stream  surface  in  this  case  suggests  that  some  dissipative 
process  ahead  of  the  body  will  occur. 
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Appendix  I 

The  main  difficulty  in  the  numerical  solution  of  the  swirling  pipe  flow  vortex  break¬ 
down  problem  lies  in  the  specification  of  the  boundary  conditions.  The  boundary  condi¬ 
tion  on  the  axis  of  symmetry  is  obvious  but  the  other  three  boundary  conditions  require 
a  higher  level  of  idealization. 

The  upstream  boundary  condition  at  z  =  0  is  the  one  of  main  concern.  Due 
to  the  limited  number  of  available  grid  nodes  (a  101  x  501  grid  has  been  an  upper 
limit  in  this  calculation)  the  swirl  generator,  i.e.  guide  vanes,  and  the  pipe  inlet  could 
not  be  simulated  at  the  same  time  as  the  flow  inside  the  pipe  ‘test-section’.  Hence, 
the  upstream  boundary  must  be  placed  at  some  location  downstream  of  the  inlet  and 
upstream  of  the  ‘test-section’.  This  requires  an  assumption  that  the  flow  at  the  upstream 
boundary  is  locally  cylindrical  and  independent  of  time  throughout  the  evolution  of 
the  calculated  flow.  Over  the  past  decade,  there  have  been  numerous  Laser  Doppler 
Anemometer  measurements  of  the  azimuthal  and  axial  components  of  the  velocity  (eg. 
Faler  Sc  Leibovich  (1978),  Escudier,  Bornstein  Sc  Zehnder  (1980),  Uchida,  Nakamura 
Sc  Ohsawa  (1985))  and  it  has  been  found  that  the  expressions  (16)  and  (17)  provide  a 
good  fit  to  the  experimental  distributions  well  upstream  of  the  location  where  vortex 
breakdown  occurs.  For  those  numerical  cases  where  vortex  breakdown  first  appears  well 
downstream  of  the  inlet  it  would  seem  therefore  that  these  expressions,  used  as  boundary 
conditions,  are  a  reasonable  approximation. 

In  those  cases  where  the  breakdown  bubble  migrates  upstream  a  change  in  the 
velocities  near  the  upstream  boundary  as  a  result  of  this  migration  is  a  clear  indication 
that  the  numerical  flow  has  become  'unphysical’.  In  many  cases  there  is  a  substantial 
evolution  of  the  ‘breakdown’  before  this  occurs.  The  weakness  in  this  boundary  condition 
requiring  further  evaluation,  however,  is  that  the  viscosity  changes  discontinuously  at 
the  boundary  and  correspondingly  gradients  in  v  (and  therefore  u)  change  relatively 
rapidly  at  low  Reynolds  number,  near  the  boundary. 

Of  lesser  importance,  but  with  its  own  problems,  is  the  specification  of  the  radial 
and  downstream  boundary  conditions.  At  the  downstream  boundary,  the  usual  uniform 
outflow  condition  is  specified.  In  the  numerical  experiments  performed,  so  long  as  no 
rotational  flow  enters  into  the  computational  domain  through  this  boundary,  changes  in 
the  location  of  this  boundary  demonstrated  that  the  effect  of  the  downstream  boundary 
was  negligible  in  these  calculations,  in  the  development  of  the  initial  breakdown.  How- 
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ever,  under  some  conditions,  it  was  found  that  multiple  recirculation  bubbles  formed 
(much  like  those  reported  by  Sarpkaya  (1971))  and  when  a  bubble  forms  very  near  the 
downstream  boundary  (which  often  happens  when  the  initial  azimuthal  component  of 
vorticity  is  large  and  positive  as  a  result  of  a  strong  axial  jet),  then  the  assumption  of 
uniform  outflow  is  clearly  violated  and  the  calculation  becomes  unphysical.  A  crude 
remedy  used  was  to  extend  the  computational  domain  further  downstream. 

The  radial  boundary  condition  is  a  somewhat  less  difficult  problem.  Typically, 
it  was  found  that  if  the  radial  boundary  is  placed  at  a  radius  of  at  least  two  core 
diameters,  then  the  numerical  evolution  is  weakly  affected  by  this  boundary.  A  similar 
weak  influence  of  the  radial  boundary  was  also  found  in  the  experiments  of  Harvey 
(1962).  In  the  model,  the  computational  domain  is  stretched  in  the  radial  direction  so 
that  there  is  a  concentration  of  grid  nodes  near  r  =  0  thus  allowing  proper  resolution  of 
the  dynamics  in  the  core  region.  The  radial  boundary  condition  is  taken  to  simulate  a 
limiting  streamline,  in  the  irrotational  region  of  the  flow,  which  remains  parallel  to  the 
pipe  wall.  This  is  physically  reasonable  so  long  as  the  diffusion  of  vorticity  is  not  large 
enough  to  invalidate  the  irrotational  approximation  and  the  boundary  layer  on  the  pipe 
wall  (which  is  not  computed)  remains  thin  and  attached. 

An  interesting  and  relevant  observation  to  be  made  is  the  difference  between  the 
‘downstream’  boundary  condition  in  the  confined  swirling  flow  case  and  that  in  other 
flows  in  which  vortex  breakdown  occurs.  In  these  other  cases  (e.g.  vortical  flow  over 
delta  winged  aircraft,  swirling  pipe  flows,  Ward-type  cyclone  chambers,  tornados  and 
waterspouts)  the  downstream  boundary  condition  on  the  vortex  core  is  ‘free’,  whereas  in 
the  confined  swirling  flow  the  rotating  endwall  forms  the  downstream  boundary  which 
results  in  the  physical  boundary  condition  that  v  =  0  at  r  —  0.  The  vortex  core  is 
therefore  fixed  to  the  center  of  the  downstream  boundary.  Having  thiB  fixed  downstream 
boundary  condition  on  the  vortex  core  accounts  perhaps  for  the  fact  that  the  ‘spiral’ 
form  of  vortex  breakdown  which  is  evident  in  the  other  examples  of  vortex  breakdown  is 
not  evident  at  the  values  of  HR2 /u  and  H/R  reported  in  the  confined  flow  experiments. 
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TABLE  1 :  ESTIMATES  OF  THE  WAVELENGTHS  OF  THE  OSCILLATIONS  IN  THE 
DISPLACEMENT  OF  THE  STREAMSURF ACE  DETERMINED  (i)  FROM  THE  EXPRESSION 
GIVEN  BY  (15)  AND  (ii)  FROM  THE  PERIODICITY  IN  THE  CONTOURS  OF  v2/r-  Vpdp/dr 
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FIGURE  1 :  CONTOURS  OF  (a)  y,  (b)  TAND  (c)  H  IN  THE  MERIDIONAL  PLANE  FOR  H/R  = 
2.5  AND  Re  =  1994.  THE  CONTOUR  LEVELS  ARE  NON-UNIFORMLY  SPACED,  WITH  20 
POSITIVE  AND  20  NEGATIVE  LEVELS  DETERMINED  BY  c-LEVEL(/)  =  MA X(variable)  + 
{i/20),nd  AND  c-LEVEL(/)  =  Mm  (variable)  x  (i/20),nd  RESPECTIVELY,  WHERE  md=  3  IN  (a) 
&  (b)  AND  ind  =  2  IN  (c).  THE  CONTOURS  ARE  PLOTTED  AT  f  =  1000,  BY  WHICH  TIME 
A  STEADY  STATE  HAS  BEEN  REACHED. 
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FIGURE  1c 
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FIGURE  2:  CO-ORDINATE  SYSTEM  USED  AND  A  TYPICAL  STREAM  SURFACE  v  = 

FOR  WHICH  r  =  a(z)  AND  THE  HELIX  ANGLES  FOR  THE  VELOCITY  AND  VORTICITY  AT 
z0  ARE  ao  AND  / ?0 . 


FIGURE  2:  CO-ORDINATE  SYSTEM  USED  AND  A  TYPICAL  STREAM  SURFACE  y/  -  ^ 
FOR  WHICH  r  =  ctU)  AND  THE  HELIX  ANGLES  FOR  THE  VELOCITY  AND  VORTICITY  AT 
z0  ARE  Go  AND  f$0. 
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FIGURE  4:  CONTOUR  PLOT  OF  v2/r-  '\/pdp/dr\N  THE  MERIDIONAL  PLANE  FOR  H/R  = 
2.5  AND  Re  =1994.  THE  CONTOUR  LEVELS  ARE  NON-UNIFORMLY  SPACED,  WITH  20 
POSITIVE  AND  20  NEGATIVE  LEVELS  DETERMINED  BY  C-LEVEL(/)  =  MA X(variable) 
x  (i/20)3  AND  C-LEVEL(z)  =  MIN(i /ariable)  x  (//20)3  RESPECTIVELY.  THE  CONTOURS 
ARE  PLOTTED  AT  f  =  1000,  BY  WHICH  TIME  A  STEADY  STATE  HAS  BEEN  REACHED 
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FIGURE  4:  CONTOUR  PLOT  OF  v*/r-  '\/pdp/drm  THE  MERIDIONAL  PLANE  FOR  H/R  = 
2.5  AND  Re  =1994.  THE  CONTOUR  LEVELS  ARE  NON-UNIFORMLY  SPACED,  WITH  20 
POSITIVE  AND  20  NEGATIVE  LEVELS  DETERMINED  BY  C-LEVEL{/)  =  MA X(variable) 
x  (i/20)3  AND  C-LEVEL(i)  =  MIN  (variable)  x  (i/  20)3  RESPECTIVELY.  THE  CONTOURS 
ARE  PLOTTED  AT  f  =  1000,  BY  WHICH  TIME  A  STEADY  STATE  HAS  BEEN  REACHED. 
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FIGURE  6:  THE  EVOLUTION  OF  THE  AXISYMMETRIC  VORTEX  BREAKDOWN 
PHENOMENON  IN  SWIRLING  PIPE  FLOW  FOR  THE  CASE  OF  l/c=  1.5,  l/Vc  =  1.25  AND 
Re=  W ,  a/v  =  250.  THE  COMPUTATIONAL  DOMAIN  HAS  A  RADIAL  EXTENT  OF  5a 
AND  AN  AXIAL  EXTENT  OF  35a  IN  WHICH  351  GRID  POINTS  ARE  UNIFORMLY 
DISTRIBUTED  IN  THE  AXIAL  DIRECTION  AND  51  GRID  POINTS  ARE  NON-UNIFORMLY 
DISTRIBUTED  IN  THE  RADIAL  DIRECTION  AND  ARE  CONCENTRATED  NEAR  r=  0. 
THE  TIME  INCREMENT  FOR  THE  EVOLUTION  IS  <U=  0.01  AND  THE 
NON-DIMENSIONAL  TIMES  AT  WHICH  THE  CONTOURS  OF  (i)  v.  <n)  ;/  AND  fin)  HAVE 
BEEN  PLOTTED,  ARE  INDICATED  IN  EACH  FIGURE. 
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FIGURE  6  (iv)  t=80 
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FIGURE  6  (vi)  t=240 
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FIGURE  6  (vii)  t  =  280 
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FIGURE  6  (x)  t=416 
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FIGURE  6  (ix)  t  =  376 


FIGURE  7:  CALCULATED  CONTOURS  OF  [\)y>.  (ii);/  AND  (iii)»  FOR  A  FLOW  WITH  Ve  = 
1.75,  Wc=  1.6  AND  Re  INITIALLY  300.  (THE  COMPUTATIONAL  DOMAIN  IS  OF  RADIAL 
EXTENT  4a  AND  AXIAL  EXTENT  36a  WITH  361  UNIFORMLY  DISTRIBUTED  GRID 
NODES  IN  THE  AXIAL  DIRECTION  AND  41  NON-UNIFORMLY  DISTRIBUTED  GRID 
NODES  IN  THE  RADIAL  DIRECTION.  THE  TIME  INCREMENT  USED  IS  St  =  0.01 .)  (a)  t  = 
227  and  Re  =  300;  (b)  f  =  250  AND  Re  =  300;  (c)  t  =  250  FOLLOWING  A  REDUCTION 
AT  f  =  227  IN  Re  FROM  300  TO  150;  (d)  t=  250  FOLLOWING  AN  INCREASE  AT  f  =  227 
IN  THE  Re  FROM  300  TO  600. 


FIGURE  7  (a)  Re=300,  t=227 


FIGURE  8:  ESCUDIER’S  (1986)  OBSERVATION  OF  THE  FORMATION  AND 
PROPAGATION  OF  AN  AXISYMMETRIC  (BUBBLE  TYPE)  VORTEX  BREAKDOWN 
DOWNSTREAM  OF  A  GU1DEVANE  SYSTEM,  (VANE  ANGLE  <t>  =  70*  AND  THE 
REYNOLDS  NUMBER,  BASED  ON  THE  PIPE  DIAMETER  AND  MEAN  AXIAL  VELOCITY, 
IS  960.) 
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